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sidered. By means of numerical calculations we demonstrate that in the large N 
ON . limit the parquet approximation gives an excellent agreement with exact results. 
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A parquet approximation (generalized ladder diagrams) in matrix models is con- 



Dh. 1 Introduction 

CD . 

43 ■ 

^ ■ The parquet approximation (or generalized ladder approximation) was introduced by 
Landau, Abrikosov and Khalatnikov in their famous consideration of the high energy 
behavior in quantum electrodynamics JTJ. Later on this approximation has been used for 
various models of quantum field theory |2j. The results of the parquet approximation 
for weak coupling are in agreement with the renormalization group approach ||. Let us 
recall that the original aim of Landau at all was to develop a non-perturbative method 
in quantum field theory. The parquet approximation leads to a closed system of integral 
equations which have meaning not only for small but also for the large coupling constant. 
In this sense the method of the parquet approximation is a non-perturbative one. There 
are arguments [ffl, 3 that this approximation is able to catch the asymptotic behavior in 
quantum field theory at small distances even in the strong coupling regime. However, 
there is a serious problem here, see 0, because it is very difficult to estimate an error of 
this approximation in realistic models for strong coupling, since we have few information 
about the true behaviour of the model in this regime. One can control the strong coupling 
regime only for some simple models such as vector models [|J and low dimensional matrix 
models Hj in the large N limit. 

The aim of this letter is to investigate the parquet approximation for matrix models 
in the large N limit. By using numerical calculations we demonstrate that the parquet 
approximation proves to be a surprisingly good approximation for N x N matrix models 

*Steklov Mathematical Institute, Vavilov 42, GSP-1, 117966, Moscow, e-mail: arefeva@class.mian.su 
^Steklov Mathematical Institute, Vavilov 42, GSP-1, 117966, Moscow, e-mail: zubarev@class.mian.su 



1 



in the large N limit. To this end we compare the numbers of parquet diagrams with 
suitable weights with the total number of diagrams with the same weights in the large N 
limit. As the weight of diagrams we take the corresponding powers of N~ l and coupling 
constants. It is occurred that in the large N limit these two numbers that are functions 
on coupling constants and the number of external lines almost coincide for all values of 
the coupling constants including the strong coupling limit. 

The behavior of the weighted number of all diagrams in the large N limit is well known 
[0. It is given by the large N limit of the corresponding zero-dimensional matrix model. 
We find the weighted numbers of parquet diagrams for large N matrix theories with 
the quartic and cubic interactions. We find these numbers as solutions of closed sets of 
equations for parquet correlation functions in the leading order of iV -1 expansion. These 
equations are simpler than the parquet equations in |2| since in our equations w-channel 
terms do not contribute. 

As it is known the large N limit in QCD enables us to understand qualitatively certain 
striking phenomenological features of strong interactions 0-0. To perform analytical 
investigations one needs to compute the sum of all planar diagrams. At present this 
problem has been solved only in a few simple models such as zero- and one-dimensional 
matrix models and two-dimensional QCD. It was suggested || that in the large N limit 
the dynamics is dominated by a master field. The problem of finding of the master field 
has been discussed in many works []i~0"|-||22||. Recently the problem of construction of the 
master field has been solved in g|. lt 

was shown that the master field satisfies to standard 
equations of relativistic quantum field theory but fields are quantized according to a new 
rule. An explicit solution of these equations is a rather non-trivial problem and it seems 
reasonable to develop some approximated schemes to study the master field. 



There were several attempts of an approximated treatment of the planar theory [24 



In particular, there were approximations dealing with counting only some part of 
planar diagrams. In [27. 28] a so-called half-planar approximation to the planar theory was 
considered. The distinguished feature of this approximation is that one can write down 
a closed set of Schwinger- Dyson- like equations for a finite number of Green's functions. 
The half-planar approximation gives a good agreement with the known exact answers 



for the large N models for a rather wide range of coupling constant [^]. However, 
this approximation does not reproduce the correct strong coupling limit of exact planar 
Green's functions. A recursive scheme which could improve the half-planar approximation 



was proposed in |29fl. It turns out that this scheme is very closed to the approximation 
adopted in this paper. 

The paper is organized as follows. In Section 2 we give a general definition for the 
planar parquet approximation. In Sections 3 and 4 we compare the planar parquet and 
the planar solutions of zero-dimensional matrix models. We find a surprising coincidence 
of these solutions for all range of the coupling constant up to 0.1 percent. Moreover, we 
find that the planar parquet solutions for the quartic and cubic interactions have phase 
transition points which also exist in the planar solutions of the corresponding one-matrix 
models. The values of phase transition points for the planar parquet and the planar 
approximations are turned out to be very closed. 
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2 Planar Parquet Equations 



In this section we define the planar parquet approximation for the d-dimensional matrix 
model with the cubic and quartic interaction terms, 

S = [ d d xTr[-dM ■ dM + -m 2 M 2 ^M 3 + -^-M 4 ]. (2.1) 

J L 2 2 37V 1 / 2 AN K J 

Here M = (My(x)), i, j = 1,...,N is an Hermitian matrix field. The planar Green's 
functions are defined as N — > oo limit of invariant correlation functions of products of 
matrices 

U n ( Xl ,...,x n )= lim^j^^ J DMTr(M(x 1 )...M(x n ))exp(-S). (2.2) 

They satisfy to the planar Schwinger- Dyson equations |10], [15] 

(- A +m 2 ) Xl IL n ( Xl , x n ) - AII n+ i(a;i, £i, jc a , x„) + 5fil n+2 (xi, x 1 , x 1 ,x 2 , ...,x n ) - 

n 

5(xi - x m )U m _ 2 (x 2 , x m _i)n n _ m (x m+ i, x n ) = 0. (2.3) 

m=2 

These equations look almost as the Schwinger-Dyson equations for a scalar theory but 
there is an essential difference in the form of the Schwinger's terms. Planar Green's func- 
tions are symmetric only under cyclic permutations of arguments. We define the planar 
parquet approximation as an approximated perturbative solution of equations (|2.3| ) that 
takes into account only a part of full series on coupling constants. This part is specified by 
a requirement that all three- and four-point vertex parts (1PI parts) are composed from 
two-particle reducible (2PR) diagrams. Note that the perturbative definition is meaning- 
ful since the perturbative solution of the planar Schwinger-Dyson equations converges at 
least for asymptotically free interactions. So, this is the case for the interaction (271) if 



we average the signs of the couplings in the suitable way as well as for SU(oo) QCD. 

More precisely, the definition of the planar parquet approximation can be done using 
a notion of a so-called skeleton expansion. The skeleton expansion contains only a subset 
of all planar diagrams. In this subset one takes only those diagrams that contain no 
propagators, three- and four- vertices insertions. Following 't Hooft [50] we call the two-, 



three- and four-point functions the "basic Green's functions". In the four-dimensional 
space-time for the cases of interaction ( |2.1| ) as well as of Yang-Mills theory just these 
vertex functions contain divergences. To get the planar parquet approximation one has 
to take the set of skeleton diagrams and replace the bare propagator and bare three- and 
four-vertices by the parquet propagator and the corresponding vertex functions. 

The parquet basic Green's functions are defined as solutions of the set of equations that 
is presented below. The propagator is defined as a solution of the exact Schwinger-Dyson 
equation 



-<» - (2.4) 



(all contributions of tadpole diagrams are dropped out). Here 
= D(p), 
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are the full (within the parquet approximation) and the bare propagators, respectively, 

= r 3 (p,g), = T 4 (p,q,k) 

p 

are three- and four-point vertices. The three-point parquet vertex is defined as a solution 
of the following equation 



+ 



-p-q 




+ 



l^p 



<1 




(2.5) 



Equation fl2.5|) means that any 1PI 3-point diagram (except the bare one) can be decom- 
posed into three- and four-point subdiagrams connected by two internal lines so that the 
four point subdiagram is 2PI with respect to those two lines that connect the four- and 
three-point subdiagrams. It is evident that such four-point subdiagram can be either the 
bare vertex either a "dumb-bells" diagrams or a V- vertex. Here the filled vertical box V 
(horizontal box H) 



H(p,q,k) 



t 



V{p,q,k) 



is the part of the four-point vertex function that is 2PR in the t-channel (s-channel) and 
is not 2PR in the s-channel (t-channel). The two particle reducibility in the s-channel 
means that the if-vertex can be represented in the following way 



p —p—q—k 



p+q-l 



+ 




+ 



P+9- 




+ 



p+q-l 




+ 



p+q-l 



(2.6) 



Equation ( |2.6p describes the structure of if-type diagrams. It means that any H- 
diagram can be decomposed into two four-point subdiagrams connected by two internal 
lines so that at least one of these subdiagrams is 2PI in the s-channel. The similar relation 
holds for V- vertex. The vertices V and H are related by the cyclic permutation of external 
points 

V(p,q,k) = H(q,k,-p-q-k). (2.7) 
The parquet approximation for the four-vertex is defined by 



• 

—p—q—k 



+ 



(21 



Equation ( |2.8| ) states that 1PI contributions into the 4-point Green's function T 4 are 
presented by the sum of the bare 4-point vertex, horizontal and vertical vertex functions 
and diagrams that are 2PR in both s- and t-channels. 
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Equations fl2.5|)-(|2~8|) together with the Schwinger- Dyson equation for propagator ( |2.4[ ) 
form the closed system of non- linear integral equations for the functions D(p), Ts{p,q), 
r 4 (p, q, k), H(p, q, k) and V(p, q, k). Just these equation we call the planar parquet equa- 
tions for the basic functions. Note that the planar parquet equations are different from the 
usual parquet equations since in the planar approximation so-called w-channel terms 
do not contribute. 



3 Weighted Number of Planar Parquet Diagrams for 
Quart ic Interaction 

In this section we examine the weighed number of planar parquet diagrams in the matrix 
theory with quartic interaction. This means that the dependence on momenta variables 
in (|2.4|) - (|278|) must be neglected. In this case equations (|2~4])- (|2.8| ) are reduced to the 
following system of algebraic equations 

D = 1 - 2gD 2 - gD 4 T 4 . (3.9) 

T 4 = -g + H + V, (3.10) 

H = -gD 2 T 4 + VD 2 T 4l (3.11) 

V = -gD 2 T A + HD 2 T A . (3.12) 
The graphical representation of equations ( |3.9| )-( |3TT2| ) is given on Figure 1. 



XT 



+ 



x + 



* * t 



+ 




Figure 1. Graphical representation of equations ( |3.9| )- (|3~TT| ). 



Four equations (|3.9|) - (|3.12|) contain four unknown functions D(g), T(g), H(g), V(g) 
and can be solved explicitly. Excluding T 4 , H and V from (|3.9|)-( [3~T2] ) we get 



g s D b + g 2 D 5 + hg 2 D A + hgD A + (1 - 5g)D 2 -2D +1 = 0. 



(3.13) 



Below on the second line of Tables 1 we present numerical solutions of ( 3.13| ) for different 
values of the coupling constant g. Corresponding T 4 is presented on the fourth line of this 
Table. 
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Table 1. 



9 


0.001 




0.01 




0.1 




1 


10 


100 


1000 


D 


0.9980 




0.9808 




0.8576 




0.5161 


0.2129 


0.07372 


0.02400 


£)(P0 


0.9980 




0.9808 




0.8576 




0.5161 


0.2130 


0.07374 


0.02401 


r 4 


-9.980 


io- 4 


-9.813 


io- 3 


-0.087^ 


16 


-0.6896 


-5.823 


-54.38 


-531.3 


r (pQ 

1 4 


-9.980 


IO" 4 


-9.813 


10" 3 


-0.087^ 


*6 


-0.6900 


-5.835 


-54.55 


-533.1 



It is instructive to compare these results with the exact weighted numbers of the planar 
diagrams. These numbers are solutions of the corresponding planar theory. The zero- 
dimensional planar Green's functions II^ for the quartic interaction satisfy the Schwinger- 
Dyson equations 

n 

U n = -gU n+2 + n m _ 2 n„_ m , (3.14) 

m=2 

that are zero- dimensional analogs of (|2.3| ). These equations admit the following solution 

U, = g-\1-D^), (3.15) 
where the planar propagator D^ 1 ' = n 2 satisfies the equation 

27g 2 D ipl) 2 + (1 + 18g)D w - 1 - 16g = 0. (3.16) 

The numerical values for the planar propagator and the four-point vertex = 
( jD (p0)-4p 4 _ 2(L>M)2] are a i so presented on Table 1. 

We see that numerical values of Green's functions in the planar parquet approximation 
and in the planar approximation for equal values of the coupling constant are strikingly 
closed, namely the difference is less then 0.1 percent for all range of g. 

It is also instructive to compare an asymptotic behaviour of Green's functions in the 
planar parquet approximation with the planar approximation. From equation (|3.13|) we 
have the following asymptotic behavior of D(g) in the strong coupling regime 

D(g) ~ a var g- x '\ (3.17) 

where a par is the solution of an equation 

a lar + 5 «par ~ 5 «Jar + 1 = °" ( 3 - 18 ) 

The numerical solution of equation ( 3.18Q is 

a par = 0.76948 . (3.19) 
The asymptotic behavior of the planar propagator follows from equation ( |3.16| ) 

D^ l \g) ~ a pl g- l '\ (3.20) 

where 

a pl = = 0.76980 . (3.21) 
9 

The agreement of ( 3.19 ) and ( 3.21| ) is more then susceptible. 
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We also compare propagator and the four- vertex function in the parquet approximation 
with the exact results for negative values of the coupling constant (see Table 2). 



Table 2. 



9 


-0.001 


-0.01 


-0.05 


-0.08 


-0.0833 


-0.0834 -0.0864 -0.0865 


D 


1.0020 


1.0210 


1.1332 


1.2959 




1.3889 


£)(P0 


1.0020 


1.0210 


1.1332 


1.2963 


1.333 




r 4 


0.0010 


0.0102 


0.05806 


0.1207 




0.1728 


r ipi) 

1 4 


0.0010 


0.0102 


0.05806 


0.1214 


0.1403 





Remarkably that the solution of the parquet equation has a phase transition for g = g , 
—0.0865 < go < —0.0864 that is related with the root branch point of the solution. In 
the planar approximation the phase transition point is caused by the square root branch 
point and it is equal to go = —1/12 = —0.0833(3). 

4 Weighted Numbers of Planar Parquet Diagrams 
for Cubic Interaction 

In this section we consider the weighed numbers of planar parquet diagrams for the matrix 
model with the cubic interaction. For this case planar parquet equations (|2.4|) - (|2.8|) reduce 
to the form 

D = 1 + \D 3 T 3 , (4.22) 
T 3 = A + D 2 VT 4 + D 3 T 3 3 (4.23) 
r 4 = H + V + D*Tj, (4.24) 
H = D 2 T 4 V + D 3 T 2 3 V + D 3 T 4 Tl (4.25) 
V = D 2 T A H + D 3 YlH + D 3 T A Tl. (4.26) 
A graphical representation of these equations is given on Figure 2. 
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+ 



+ 








Figure 2. Graphical representation of equations fl4.22|) -( 4.25|) . 
Note that using equations ( |4.24 )-( 4.26| ) one can rewrite the equation for r 3 in the form 

r 3 = A + XD 2 T 4 + XD 3 Tl (4.27) 

or graphically 

Equation (|4.27|) is just the Schwinger- Dyson equation for the 1PI three-point Green's 
function T 3 . 

Excluding T 3 , T 4 , H and V from equations (t4 . 24|) - . 2 6|) we get the equation for D 




2X e D 9 - 3\ 4 D«(D - 1) + X 2 D*(D - l) 2 - (D - 1 



0. 



(4.28) 



Let us compare the solution of the planar parquet equations (|4.24j) - (|4.23|) with the 
solution of the planar theory. The planar 2-point and 1PI 3-point Green's function are 
given by formulas ||: 



l-3r 



;i - 2r) 

where r is the subject of the relation 



2' L 3 ~ A 



;i-2r) 2 (l-4r) 
fl-3r) 3 



r(l-2r) 2 . 



(4.29) 



(4.30) 



The results of numerical solutions of planar parquet equations ( |4.24| )-( |4~23| ) as well as 
equations ( |4.29|) -( |4.30|) for planar Green's functions are presented on Tables 3 and 4. 



Table 3. Perturbative solution. 
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A 


0.001 


0.01 


0.1 


0.2 


0.27 


0.272 


0.273 0.279 0.280 


D 


1.00 


1.00 


1.0104 


1.0486 


1.1197 




1.1458 




1.00 


1.00 


1.0104 


1.0486 


1.1200 


1.1245 




r 3 


0.0010 


0.010 


0.10106 


0.21084 


0.31577 




0.34720 


f (pO 

1 3 


0.0010 


0.010 


0.10106 


0.21084 


0.31644 


0.32202 





Table 4. Nonperturbative solution. 



A 


0.001 


0.01 


0.1 


1 


10 




100 


1000 


D 


-26783 


-864 


-30.78 


-1.96 


-0.2? 


36 


-0.0565 


-0.0119 


£)(P») 


-2.51 • 10 5 


-2640 


-39.6 


-2.00 


-0.2? 


37 


-0.0566 


-0.0120 


r 3 


1.39 • 10- 6 


1.34- 10~ 4 


0.0109 


0.0.391 


5.48 




58.6 


594 


r (pQ 

1 3 


1.58 ■ 10~ 8 


1.43 ■ 10~ 5 


0.00653 


0.375 


5.42 




58.2 


590 



Table 3 corresponds to a perturbative solution. When A = we have D — 1 and 
T 3 = that is related to the free theory. There is a phase transition point A = Ao 
where 0.279 < Ao < 0.280 for the planar parquet approximation and Ao = v^2/(3a/3) = 
0.2722 for the planar approximation. Table 4 corresponds to another branch that is a 
nonperturbative solution. 

From (|4.28|) we find the following asymptotic behavior for D in the strong coupling 
regime 

D(g) ~ c par A~ 2/3 , (4.31) 

where c par is satisfied the equation 

2c 9 par + 3c 6 par + 1 = 0. (4.32) 



The numerical solution of equation Q4.32 ) is 



c par = -1.18823 . (4.33) 



An asymptotic behavior of the planar propagator is found from ( (4.29| ), ( |4.30| ) and has the 
form: 

D( pl \\) ~ c pl \- 2/3 (4.34) 

with 

3(4W 3 

c p i = — — = -1.19055 . (4.35) 



5 Conclusion 

In conclusion, in this paper we found that the planar parquet approximation gives a 
unexpectedly good agreement with the known exact results for large N matrix models. 
The parquet approximation in quantum field theory is one of few known approximations 
which admits a non-perturbative closed set of equations for finite Green's functions. 

The results of this paper mean that the class of parquet diagrams is rather com- 
prehensive among of all planar diagrams of matrix theories. This makes reasonable an 
investigation of planar parquet approximation for QCD that will be the subject of our 
further work. 
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